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Conformations, hydrodynamic interactions, and instabilities of sedimenting semiflexible filaments 
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The conformations and dynamics of semiflexible filaments subject to a homogeneous external (gravitational) 
field, e.g., in a centrifuge, are studied numerically and analytically. The competition between hydrodynamic drag 
and bending elasticity generates new shapes and dynamical features. We show that the shape of a semiflexible 
filament undergoes instabilities as the external field increases. We identify two transitions that correspond to the 
excitation of higher bending modes. In particular, for strong fields the filament stabilizes in a non-planar shape, 
resulting in a sideways drift or in helical trajectories. For two interacting filaments, we find the same transitions, 
with the important consequence that the new non-planar shapes have an effective hydrodynamic repulsion, in 
contrast to the planar shapes which attract themselves even when their osculating planes are rotated with respect 
to each other. For the case of planar filaments, we show analytically and numerically that the relative velocity is 
not necessarily due to a different drag of the individual filaments, but to the hydrodynamic interactions induced 
by their shape asymmetry. 


I. INTRODUCTION 

Semiflexible filaments are fundamental constituents of 
micro-biological systems, where microtubules and actin fila¬ 
ments serve as scaffolds for cellular structures and as routes to 
sustain and guide cellular transport systems. HI Microtubules 
are also the main structural elements of cilia and sperm flag¬ 
ella, where their relative displacement and deformation due 
to motor proteins gives rise to the flagellar beat and hydrody¬ 
namic propulsion. Cl El Microtubules and flagella can be seen 
as elastic filaments interacting with their own flow field. The 
ability to visualize, assemble, and manipulate biological and 
artificial semiflexible polymers ffia poses new fundamental 
questions on the dynamics of filaments when elastic and hy¬ 
drodynamic forces compete. 

The dragging of stiff rods through a viscous fluid has been 
studied in detail. JS) A single rod does not reorient, but falls 
with its initial orientation. A more complex dynamical be¬ 
havior can be expected and is indeed observed for semi¬ 
flexible filaments when the curvature or stretching elasticity 
competes with the hydrodynamic interactions. EHm Single 
dragged semiflexible filaments bend into a shallow V-shape 
to balance the higher drag at both ends ca and their end- 
to-end vector aligns perpendicularly to the external held. IfTTl 
For strong drag, higher modes have been found to be excited; 
this generates VT-shapes initially, which then relax back into 
horseshoe-like tZ-shapes. na Here, the dynamics seems to be 
constrained to the plane initially defined by the direction of 
the external held and the filament itself. However, these in¬ 
vestigations address the problem from a deterministic point of 
view, and little attention has been paid to the dynamic stability 
of the resulting shapes. In all cases, the dragged and deformed 
semiflexible filament initially defines the settling plane, but 
the stability of the filament’s planar shape has not been inves¬ 
tigated as function of the external held or the relative position 
of possible neighboring filaments. 
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Here, we focus on the full three-dimensional shape of one, 
two, and three semiflexible filaments sedimenting in a ho¬ 
mogeneous external held. We incorporate the hydrodynam¬ 
ics into the equations of motion for the filament shape via 
the Oseen tensor, valid in the limit of zero Reynolds num¬ 
ber. As a result of our numerical and analytical analysis, we 
And that the deformations confined to a plane become unsta¬ 
ble with respect to normal perturbations at a threshold value 
Z?j of the strength B of the external held, which is smaller 
than the threshold B 2 where initial, transient VT-shapes be¬ 
come excited, see Fig. [T] Thus, with increasing strength of 
the external held, two instabilities and transitions to new sed¬ 
imentation modes are predicted. The first transition is from 
a stable planar (/-shape with little bending to a stationary 
horseshoe-like (/-shape with out-of-plane bending. The sec¬ 
ond transition at stronger fields excites a metastable W-shape, 
also with out-of-plane bending, which then “relaxes” into a 
non-stationary asymmetric (/-shape. As result, there exist two 
families of shapes, where the elastic forces are balanced by a 
conformation-dependent drag. 

We consider next the interaction between two filaments 
in an external held. Indeed, while the dynamics of an iso¬ 
lated filament is an indispensable knowledge needed to under¬ 
stand the case of n > 1 interacting filaments, many situations 
are characterized by elastic slender objects interacting via the 
generated flow field; cilia, ||71 [13] sperm, lfT4lfT5l and E. Coll 
bundles US] El are probably the most relevant from a biolog¬ 
ical point of view. 

It is known that the sedimentation behavior of colloids 
can be quite complex. The interaction of sedimenting par¬ 
ticles has been studied in considerable detail for spherical 
colloids. lITSl [ 19 ] Two particles sediment together, but don’t 
follow the direction of the external held, and move instead 
under an angle with respect to it. For more particles, many 
different dynamical behaviors can be found, in particular pe¬ 
riodic motions where particles “dance” around each other. lfT9l 

For dragged semiflexible filaments, the dynamical behavior 
is even more complex. iHOl In particular, we show that two fil¬ 
aments (Fig. [T]| attract each other, repel each other, or spin 
around the held depending on the intensity of the external 
held. 
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Figure 1: Snapshots from simulations of single filaments dragged by the external homogeneous field B = mgL^/K, where L is the filament 
length, g the external field, and K the bending stiffness. Left: For weak field (B < fij) the filament bends into a F-shape (in dots), dominated 
by the X2,z mode. Center: As the field strength increases, higher modes with an out-of-plane component are excited, and the filament drifts 
sideways. Right: For even stronger fields {B > B 2 ) further symmetries are spontaneously broken, and the filament rotates following a helical 
trajectory. Corresponding movie is shown in the ESL The vertical distance between the frames is reduced to enhance the visualization of the 
filament conformations. 


We focus here on the stability of the sedimentation plane 
for different field intensities and on the origin of the relative 
velocity. In particular, we want to see whether the velocity 
difference is due to different shapes or to the broken up-down 
symmetry. For even more filaments, the dynamics becomes 
unsteady at much weaker external field strength than expected 
from the two-filaments case. 


II. MODEL AND METHODS 
A. Discrete Model 

In the simulation, filaments of length L = b{N — 1) are rep¬ 
resented by N mass points of mass m connected by harmonic 
bonds of length b, with the potential 

= y E\|Rd-^)^ (1) 

^ !=0 

where R, = r,+i — r, is the vector connecting the consecutive 
points, i G {0,... ,N}, and kb is the force constant. To account 
for filament stiffness, we introduce the bending potential 


a particle 

Fp = -mge„ (3) 

with the unit vector e, along the z axis of the Cartesian coor¬ 
dinate system (cf. Fig.[^. 

Inter- and intrafilament hydrodynamic interactions substan¬ 
tially influence the filament dynamics. Hence, we apply the 
equation of motion 

7 o(r,-u,)=Ff (4) 

for the overdamped dynamics. 11211 where, u, is the back¬ 
ground flow velocity at the position r, of the particle, 70 = 
3nrib is its friction coefficient for a fluid of viscosity 77 , and 
Fp is the sum of all conservative forces. ifTOl We compute the 
background flow velocity explicitly via the Oseen tensor. EH 
Thus, the final equations of motion are 

r, = EH,;FC, (5) 

J=o 

with the hydrodynamic tensor 


U = ( 2 ) 

with the bending rigidity k.lj^O) In addition, the mass points 
are exposed to an external gravitational held with the force on 


H ,7 = H(r,--y) 
= (1-5,;) 


3b 




.^Yonj 

and nj = n - Tj, nj =1 Tij I, and f,y = nj/nj. 


+il 

7o 


( 6 ) 
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If needed, excluded-volume interactions are implemented 
via a short-range Yukawa potential between points of differ¬ 
ent filaments, which implies a minimal effective distance dur¬ 
ing the simulations. A low-amplitude white noise is added 
to avoid metastable states. The noise is not considered to be 
of thermal origin as it is chosen to be negligible compared 
to the other hydrodynamic and mechanical forces and barely 
influences the stationary settling dynamics for the considered 
external fields. 


B. Parameters/Methods 

We set the hydrodynamic diameter of a point equal to the 
bond length, as in the Shish-Kabab model of Refs. ifTOl [TJ] 
ED, thereby fixing the aspect ratio to b jL. Lengths are mea¬ 
sured in units of the hydrodynamic diameter b and time in 
units of job^ jx. This choice eliminates the friction coeffi¬ 
cient and the bending rigidity from the equations of motion. 
The force constant for the bonds is set to k^b^ /k — 1, result¬ 
ing in a maximum extension of ±0.6% of the length over the 
investigated range of parameters. In these units, the external 
field strength mg becomes G = mgb^ /k. For convenience and 
an easier comparison with results of Ref. Qol , we characterize 
the external force hy B — N^G, ox B — mgl} jx. For B ^ 1, 
the bending rigidity dominates and the filament is essentially 
straight. We consider only filaments of length L = 301? in the 
following. For excluded volume interactions, the minimal ef¬ 
fective distance is approximately 5b. The equations of motion 
are integrated with an adaptive time-stepping Velocity-Verlet 
algorithm. Il22ll^ 


C. Continuum Model 


For an analytical description of the filament dynamics, we 
adopt a continuum model. The equation of motion of the point 
r''(s,t) (—L/2 < s < L/2) along the contour of filament v is 
given by ll^ 

rL/l 

dy{s,t) =Y, ds'H{r''{s)-r^{s'))f^{s'), (7) 

yj-L/2 


where is the external force density and the index jj. indicates 
the various filaments. As before, the hydrodynamic tensor 
H(r''(s) — r^(s')) comprises the Oseen tensor and the local 
friction. Explicitly, it reads as 


H(R) 


1 Q(|R|-^) 

STrr; |R|3 


[IR2±RR^] 


g(R) j 
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( 8 ) 


Here, 0(.r) is the Heaviside function, y— 37X7] is the friction 
per unit length, andR = r''(i)— r^(i').||^|25l The force den¬ 
sity f comprises bond, bending, and gravitational forces. In 
the limit of a rather stiff filament, it can be written as 

= l,k,T - Tj) r"(,) +fj(,), (9) 


with the persistence length L. lE^ |27]| In the following, we 
will neglect the bond term, i.e., the term with the second 
derivative and focus on bending stiffness only. 

The expansion 


n=0 

in terms of the eigenfunctions of the biharmonic operator, 
i.e.. 


IpksT = —(pn (11) 

O S Tfi 

with suitable boundary conditions. Il26l ED leads to the equa¬ 
tions of motion of the mode amplitudes 


dtX 


V 

n 


LEHrr 

M /=0 


7 

T/ 


xfit)- 




( 12 ) 


The matrix representation of the hydrodynamic tensor is 

rL/2 

H,7 = / dsds' (13) 

J-L/2 

The eigenfunctions 0„(i) and relaxation times T„ are well 
known. ||5] For convenience, we summarize them in 

Appendix [A| However, Eq. ( [T2| is nonlinear and thus cannot 
straightforwardly be solved. Eor the current analysis, the sec¬ 
ond and forth mode are most important; they are responsible 
for the V— and W-shape displayed in Eig.[^ 

To characterize the numerically obtained filament confor¬ 
mations, we calculate the mode amplitudes 

1-1/2 

Xn{t)= dsr''{s,t)p„{s). (14) 

J-L/2 

The components of the vector indicate the importance of 
the mode in the Cartesian directions. Eor example, the mode 
amplitude X 2 z measures how much the filament bends along 
the z direction into a V-like shape. 

As for the discrete model, we scale lengths by the fila¬ 
ment diameter b and time by yob^/ksTlp. The latter is iden¬ 
tical to the time scaling of the discrete model, since x = 
ksTlp. This implies for the strength of the external force 
G = pgb^/ksTlp =mgb^lx. 


III. RESULTS 

A. Deformation and Dynamics of Single Filament 

The filament is initially oriented along the x axis of the 
reference frame. After a certain time, the dragged filament 
reaches a stationary shape and velocity. Examples of confor¬ 
mational sequences for various field strengths are displayed in 
Fig.[T] We characterize the shapes via Eq. ( [T4] l in terms of the 
mode amplitudes. In Eig.[^ the most important stationary am¬ 
plitudes are presented. Below a critical field B\ ~ 1200, the 
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Figure 2: Stationary mode amplitudes of a single semiflexible fila¬ 
ment as function of the external field B. The shaded areas indicate 
the 66% confidence interval. When B < only planar modes are 
excited, and the filament stays in the plane defined by its initial orien¬ 
tation and the orientation of the applied field, here the xz plane. For 
B > Bj, an out-of-plane mode Xl.L is excited. For B > B^, the out- 
of-plane component X2.±^ ths bending component Xl.z saturates, and 
the amplitude Xa,z becomes important (visualized in Fig.fTl. In black 
crosses indicate the maximum value of X4,z before it decays. The 
resulting shape is asymmetric and spirals around the z axis with fre¬ 
quency I £0 1 / 0)2 (orange line), with o >2 the frequency of the second 
mode (Eq. 0). 


filament shape is governed by planar modes (green and black 
lines), where Xi.z dominates and, thus, the characteristic V- 
shape appears. 

In simulations restricted to a two-dimensional plane, or 
in three-dimensional simulations without noise, ifTHl the fila¬ 
ment dynamics is localized in the xz plane and filaments bend 
into a planar W-shape for fields B > ^ 1800. In con¬ 

trast, in our three-dimensional simulations with weak noise, 
we find that the planer filament conformations are metastable 
for B*i < B < B 2 , and also modes along the y axis are excited. 
We characterize the out-of-plane filament shape and dynamics 
by the mode amplitude X2.±{t)’ where 


Z2,±(f) = X2x{t) + iX2,yif) = (15) 

In the stationary state, an (/-shaped and deck-chair-like con¬ 
formation is assumed with out-of-plane bending (see Fig. [^. 
The filament orientation is fixed and X 2 ± = X2y (blue line in 
Fig.|g. Since its shape is asymmetric, the filament drifts side¬ 
ways while settling in the external field. 

When B > B 2 , the mode X4,z becomes important at early 
times, leading to a temporary W-shape (Fig.|^. The trajectory 
for B ~ 3000, displayed in Fig. [1] shows the initial W, which 
later turns into an asymmetric (/-shape, in which one arm is 
longer than the other. The appearing shape is stable; however, 
because of its asymmetry, the mode amplitude Xi ,z is non-zero 
and the filament rotates around the z axis with frequency (O, 
see Fig. [^(orange line), which we determined via Eq. ( [T5] l. 

In Fig. we characterize the helical trajectories by the 
pitch, radius, and rotation frequency (B > Bp. As the field 
increases the rotation frequency increases, and the helix be¬ 
comes more tight because the radius decreases, and the pitch 


Figure 3: Pitch, radius, and rotation frequency of the helical tra¬ 
jectories when B > B 2 - The radius and pitch seem to diverge in the 
proximity of the transition. Bars are standard deviations. 


shortens. The ratio between the pitch and the radius defines 
the helix angle a « \%/9, constant for all B > B"^- Approach¬ 
ing the transition point from above, the rotation frequency 
vanishes, and both radius and pitch diverge because the tra¬ 
jectory straightens. 

In contrast, in the deterministic dynamics of previous 
studies, ca the VT-shape was found to decay only into the sta¬ 
ble and symmetric planar horseshoe shape. 


B. Conformations and Dynamics of Two Interacting Filaments 

1. Weak Field - Relative Velocity 


As shown in Sec. Ill A the stationary shape of a single fil¬ 
ament in weak fields B < Bj is of V-shape, which breaks the 
bottom-top symmetry. This is sufficient to generate an effec¬ 
tive attraction between sedimenting filaments with the same 
shape. To characterize this interaction, we compute the rela¬ 
tive velocity Av between the centers of mass of two filaments 
of equal shape along the sedimentation direction. The fila¬ 
ments remain localized in the xz plane and are separated by a 
distance d. As shown in Fig. the relative velocities exhibit 
a significant dependence on the filament separation. We espe¬ 
cially find that Av ~ for distances larger than the filament 
length. 

The distance dependence can be understood by the theoret¬ 
ical model introduced in Sec. Ill Cl For the considered filament 
shapes, 


r'('ST) = {XiAf)^ii^)^^^X2,zh{s)f (16) 


and A{s,t) = r* (s,f) +de^, the general expression for the ve¬ 
locity difference derived in Appendix [B] yields 


AVcm 


2 

X2,z^2 


(17) 


in the limit B ) 0 °. Evidently, the filaments attract each other 
due to the top-bottom asymmetry of their shapes. In the sim- 
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Figure 4: Simulations of two filaments with the same imposed shape, 
kept constant during the simulation (B = 195). The shapes are cre¬ 
ated with the given X2.z- The filaments lie in the same plane, parallel 
to the external field. The relative velocity Av scales as d^'^. The 
black lines correspond to the prediction of Eq. l |17^ , save for a com¬ 
mon factor 5 ~ 11. The theory describes correctly the trend on d, 
and the trend on X2,z holds up to Xl.z = 8 x 10^^. 


ulations, the filament shapes are determined initially by im¬ 
posing the amplitude X2,z’ which is then kept fixed. The sim¬ 
ulation results of Fig. are in agreement with our theoreti¬ 
cal prediction down to roughly the filament length. The 
power law is indeed a universal scaling, unaffected by the fila¬ 
ments shape and external field as evident from the theoretical 
considerations in Appendix [B] The dependence of Av on X2,z 
(Eq. (fTT])) is also verified for very small bending. 


2. Weak Field - Stability 

We now relax the imposed shape constraint and con¬ 
sider collective effects for two filaments, which are initially 
straight, oriented along the x axis, and displaced along the 
z axis by a distance d (cf. Fig. |^a)). For easier compari¬ 
son with Ref. mniiiii, we employ the dimensionless number 
A = (Aupper — Aiowei)/(T/2) to quantify the bending asym¬ 
metry, where Aupper, lower is the total z extension of the up¬ 
per/lower filament. As indicated in Fig. the filament curva¬ 
ture changes with time and the upper filament is bent stronger 
than the lower one. Figures (a),(b), show the curvature 
asymmetries A and the relative velocities for various exter¬ 
nal field strengths. A decreases with increasing distance d, 
indicating more similar shapes at larger distances. Hydrody¬ 
namic interactions lead to an attraction of the two filaments 
(''upper > ''lower), in agreement with the imposed-shape ap¬ 
proximation studies of the last section. Indeed, the constant- 
shape approximation still gives the correct {L/d)^ power-law 
dependence for d/L^ 1, while the magnitude of the defor¬ 
mation, X2z^'^’ fitted. When the filaments approach 

each other, the generated flow field depends on the details of 
their shapes that, in turn, depends on the external field, hence 
we expect a non-universal behavior. Note that in contrast to 
Ref. cni, we find that the upper filament bends more than the 


planar rotated 



Figure 5: Snapshots of two-filament conformations for B = 195, in 
time intervals At. (a) Co-planar sedimentation. Note that the upper 
filament is more bent than the lower filament, and = 0.13. 

Axes to scale, z position translated, (b) The two filaments approach 
each other after initialization in a rotated configuration. Both fila¬ 
ments spin around the z axis, with the upper filament spinning faster 
(seeFig.|^(c-d)). 


lower filament (see also Fig.[^. 

The planar configuration of a filament is also stable with re¬ 
spect to filament rotations around the field axis, see Fig.|^(b). 
Filaments that are initially displaced along the z axis (as in 
the previous case) and rotated with relative orientation angle 
0 around the external field axis spin until the relative angle 
vanishes, as illustrated by Fig. (c). Also in this case, the 
upper filament drifts and rotates faster than the lower one, see 
Fig. [^(d). The relative velocity is essentially the same as in 
the planar case. 

Thus, two filaments sedimenting in weak fields relax to¬ 
ward a stable planar configuration one behind the other. The 
shape of the filaments is dominated by the second mode, 
pointing downwards, as shown in Fig. This mode domi¬ 
nates and it breaks the mirror symmetry of the hydrodynamic 
interactions even for filaments of the same shape. Note that, in 
contrast to the single filament case, the system does not reach 
a stationary state velocity or shape, since the upper filament is 
always faster than the lower filament until the filaments touch 
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Figure 6: a)-b) Bending asymmetry A and relative velocity Av of two filaments as function of the filaments distance for Ljb = 30. The two 
filaments are in the same plane, parallel to the external field and parallel to each other. Each color corresponds to a different external field B, as 
indicated. The velocity vq is the terminal velocity given hy the resistive force theory for a rod. When d/L 3> 1, the relative velocity scales as 
d^'^. Note that filaments attract, i.e. time progresses from right to left. c)-d) Rotation angle 6 and relative velocity Av of two initially rotated 
filaments around the field axis by 0 = 18°. The relative velocity is essentially unaffected by this change. Notably, the filaments spin toward 
each other decreasing the relative angle. 


each other. 


3. Strong Field 

For strong fields, we consider two filaments, which are ini¬ 
tially displaced by 6L along the field direction. We measure 
the shape eigenvalues when the distance is 5L, in the quasi- 
stationary regime, and find that the eigenvalues exhibit the 
same behavior as those of a single filament. This means that 
for B > Bj the dynamics of each filament is dominated by 
the local flow field and not by the interactions with the other 
filament. Indeed, we find no correlations between the orienta¬ 
tions of the out-of-plane components of the two filaments for 
B > B]; the two filaments can by chance bend out-of-plane 
and drift in arbitrary directions. 

When B > B^, the filaments undergo the same transitions 
as a single filament: each of them reaches the same stationary 
shape and rotation velocity as an isolated filament. We find 
no correlations between the rotation directions of the two fila¬ 
ments: some filaments spin in the same direction, others in op¬ 
posite directions, with no preference. This highlights the rel¬ 
evance of hydrodynamic interactions between two filaments 
for external fields weaker than Bj. At stronger fields, the 
non-planar configurations generate forces that compete with 
and dominate over hydrodynamic interactions among the fil¬ 
aments and the emergent behavior is, essentially, the same as 
that of an isolated filament. 



Figure 7: Three semifiexible filaments and trajectory of one bead 
(thick line), for the external field 6 ~ 60 ^ Bj. In this case, the fil¬ 
aments form a bundle, but the relative positions change periodically. 
Inset: plot of X2,z for th® three filaments. Since they have the same 
period and constant phase shift, this is the result of a cooperative 
behavior. Corresponding movie is shown in the ESI. 


C. Three Filaments 

Given the complex dynamics of two interacting filaments, it 
is interesting to consider also the collective behavior of several 
filaments. We find in simulations of systems with more than 
two filaments an intriguing collective dynamic behavior even 
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for very weak fields (B ~ 60) and in the absence of noise. 

We focus here on the case of three filaments, see Fig.|^ For 
most (randomly chosen) initial configurations, the nearest two 
filaments form a bundle that settles faster than the third fila¬ 
ment that is then left behind. However, we find also some ini¬ 
tial configurations where all three filaments attract each other 
and form a bundle. In this case, the relative positions are not 
stationary; instead, the filaments follow a periodic trajectory, 
see Fig. (inset). In the inset of Fig. we show also that 
the shapes of the three filaments are not stationary. The mode 
amplitude X2.z of s^ch filament changes periodically, with a 
constant phase shift between them. 

Our results for one and two filaments indicate that trig¬ 
gering of a time-periodic bifurcation requires strong fields. 
However, the three-filaments results suggest that systems with 
more filaments display a very complex dynamics even for 
weak fields due to complex hydrodynamic interactions. 


IV. DISCUSSION AND CONCLUSIONS 


B < lO'*. We can estimate these parameters for biopolymers 
like actin or microtubules. Actin has a persistence length of 
Ip ~ 17/iOT, an average length L ~ lOjim, and the bending 
rigidity ff ~ 60 x The external gravitational 

field, corrected for buoyancy, is about G « 10 which 

implies Bgravity — 10^^. Microtubules, on the other hand, are 
stiffer, longer, and heavier with Ip ~ 1mm, L ~ lOO/lm <C Ip 
and G « 10~^^. ll28l This yields the effective field strength 

Bgravity — 10^*. An experimental test of our predictions is 
therefore within reach of modem centrifuges with accelera¬ 
tions of about lO^g. 

Appendix A: Eigenfunctions of a Semiflexible Filament 
The eigenvalue equation ( [TT] ) with the boundary conditions 

I I 

(‘^) |.^=± t ./2 | 5 =± Z ./2 0 ( 2k 1 ) 


We have investigated the dynamics and stability of semi¬ 
flexible filaments exposed to an external homogeneous field 
and interacting only via hydrodynamic fluid fields. Due to the 
competition between hydrodynamic interactions and bending 
stiffness, the appearing dynamical behavior is richer than for 
entropy-dominated polymers or interacting rods. 

We have shown that, for weak fields B < B*^, co-planar 
configurations of two filaments are stable upon perturbations 
that rotate the shapes relative to each other around the field 
axis. With simulations of fixed shape filaments, we have high¬ 
lighted that a V- or t/-shape is sufficient to break the hydrody¬ 
namic symmetry at low Reynolds numbers, leading to a rela¬ 
tive velocity that scales with distance as (L/d)'^. Hence, the 
difference in drag coefficients between filaments is not nec¬ 
essary to explain the faster settling velocity of the upper fila¬ 
ment. 

For external field strengths exceeding the critical value 
the hydrodynamic interactions bend the filament out of its 
principal plane. Simulations of a single bent filament show 
that the hydrodynamic forces balance the elastic force, stabi¬ 
lizing the out-of-plane shape. The resulting trajectory shows a 
drift in the direction of out-of-plane bending, superimposed to 
the settling motion. This is a novel result, not to be confused 
with the previously reported metastable W-state ca that is 
excited when B > B^. A careful analysis of the eigenmodes 
indicates that the decay of the metastable state does not, in 
general, lead to the reported planar horseshoe shape, but also 
excites an average rotation mode with respect to the field axis 
(;Ki -) and our out-of-plane bending mode X 2 ,±- The filaments 
spin then around the field axis. 

Finally, we have demonstrated that three filaments display 
an unexpected periodic dynamics even at field strengths far 
weaker than . This is in contrast to the dynamics of a pair 
of filaments that either displays a monotonic dynamics that 
relaxes the attractive force (when B is weak) or a dynamics 
dominated by the single-filament (when B > B\). 

The interesting external fields B are in the range 10^ < 


yields the eigenfunctions 

It, ( \ - ^ sinC„s \ 

" Vl \smhC„L/2 + sinUn) ’ 

J_ / cosh COS^nS \ 

Vl \coshi^„L/2 cos CnLl2j ’ 


n > 1, odd, 
(A2) 

n > 1, even. 

(A3) 


The wave numbers are approximately given by = {2n — 
l)7r/2L (n > 1), and the corresponding relaxation times 


167L4 

lpn^kBT{2n-lY' 


(A4) 


More precise eigenfunctions are provided in Ref. IMIEtI. 
The set of functions is complemented by the eigenfunction 
of the center-of-mass translation ETlI 


00 = 


1 

Vi 


(A5) 


and that of rotation of the rodlike object 


01 



(A6) 


with the relaxation time 


24-kBT' 


(A7) 


Appendix B: Relative Velocity of two Filaments 

We derive here an equation for the relative velocity between 
the centers of mass of two filaments. We restrict our analy¬ 
sis to the case of small bending amplitudes, that is equivalent 
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to consider small external fields, and filaments of identical 
shape. 

Since /f'02 (l>n{s)ds = y/Ldnfl for the exact eigenfunctions, 

the difference in the center-of-mass velocity Avcm = 
of two isolated filaments is given by 


1 

Avcm= 7 / iis<9,[r'(s,f)-r^(s,f)] 

L J-l/2 

= ^d,[xl{t)-xl{t)\. (Bl) 

Substitution of Eq. ([T^ yields 


s/LAycm= EHq,' 

n 

in 

+ f«G 

rt 

y v2 , n2 

~~Xn^\G 

n 

r„2 

Ah 

in 

O 

+ 

-lh^: 

n 

y y 1 1 f 1 

%n 


The first two terms on the right-hand side account for self¬ 
interactions of the individual filaments, the other two terms 
for the hydrodynamic interactions between the filaments. 

We simplify our considerations by assuming identical 
shapes of the filaments, i.e., we set xh — Xn ■— Xn- Moreover, 
for the constant external force the relation applies = fog^on 
independent of the particular filament. Hence, its contribution 
vanishes, which yields 


Avc,„ = 


1 

7 l 


i « 

n=\ 




(B2) 


We are primarily interested in the distance dependence of the 
relative center-of-mass velocity. Hence, we additionally ne¬ 
glect the dyadic term in the hydrodynamic tensor ([^. More¬ 
over, the local friction term vanishes in Eq. ( |B2| i, and the hy¬ 
drodynamic tensor can be written as 


HV/J _ ^ ^n{s)^o{s’) j I I 

Snri J-l /2 \r''{s) — ri^{s')\ 


(B3) 


Using the eigenfunction expansion Eq. ( [T0| , we obtain 
r''(i) -r'"(/) =Ar^^ + JJXn 

n=l 

=ArZ+Z{s,s'). (B4) 


With this definition, we obtain for AHq^ = Hq,' — 

1 

J ^^^dsds'U^)U^') (B5) 

1 1 

In the limit d = |Ar^^| 3> |E(i,i')|, Taylor expansion yields 

AHi,^ =—j^^^dsds'Us) (B6) 

and hence. 


Ay cm — 


1 1 


/■L/2 


, dsds<^n{s) 


Z{s,s')-Arl 


,21 

cm ^ /J\ 


Anr\ Vl d^ T„ " J-l /2 


W). 


(B7) 


Substituting x = sjL and setting Y—3nri Il24l . yields 


Av — 

AA T rm — 


3L2 ~ 1 fi/2 ^_^^E(x,y)-Ar2i 


4 d^ “j T„ y/T. J-l /2 


dxdx'(f)„{x)- 


(B8) 


Thus, the relative velocity decreases quadratically with the 
distance between the filaments. There is evidently no velocity 
difference when Ar^j, is perpendicular to E(s,s'). In particu¬ 
lar, there is no force between two specifically aligned rods as 
long as their director Xi is perpendicular to Ar^^. 
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